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ABSTRACT 

We present the results of Chandra observations of 13 optically-selected clus- 
ters with 0.6 < z < 1.1, discovered via the Red-sequence Cluster Survey (RCS). 
All but one arc detected at S/N> 3; though three were not observed long enough 
to support detailed analysis. Surface brightness profiles are fit to /3-models. 
Integrated spectra are extracted within R25oo; and Tx and Lx information is 
obtained. We derive gas masses and total masses within R2500 and R5oo- Cosmo- 
logically corrected scaling relations are investigated, and we find the RCS clusters 
to be consistent with self-similar scaling expectations. However, discrepancies ex- 
ist between the RCS sample and lower- 2; X-ray selected samples for relationships 
involving Lx, with the higher-z RCS clusters having lower Lx for a given Tx. In 
addition, we find that gas mass fractions within R2500 for the high-z RCS sample 
are lower than expected by a factor of ~ 2. This suggests that the central entropy 
of these high- 2; objects has been elevated by processes such as prc-hcating, merg- 
ers, and/or AGN outbursts, that their gas is still infalling, or that they contain 
comparatively more baryonic matter in the form of stars. Finally, relationships 
between red-sequence optical richness (Bgc,red) and X-ray properties are fit to the 
data. For systems with measured Tx, we find that optical richness correlates 
with both Tx and mass, having a scatter of ~ 30% with mass for both X-ray 
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and optically-selected clusters. However, we also find that X-ray luminosity is 
not well correlated with richness, and that several of our sample appear to be 
significantly X-ray faint. 

Subject headings: cosmology:observations — X-rays:galaxies:clusters — galaxies:clusters:general 



Introduction 



The Extended Medium Sensitivity Survey (iGioia et al. Ill990l . EMSS) sparked renewed 
interest in the ongoing search for clusters of galaxies at high redsh ift. Since then, nu- 



merous high redshift surveys have been carried out, both op tically (iGilbank et al.l 12004 



Donahue et al 



20021: 



Postman et al.l 



1996 



Bower et al. 



1994 



ray JValtchanov et al. II2OO4I : ISauer et al.lbooi Iwilkes et al" 



and others), and in the X- 



2000l . and others). The moti- 



vations for such searches are multifaceted, but the most compelling of these are cosmological 
in nature. 

Clusters of galaxies are an important source of information about the underlying cos- 
mology of the universe. They are considered to be essentially "closed boxes" , meaning that 
the primordial matter that they were initially assembled from has remained trapped in their 
deep potential wells since they were formed. This makes them ideal objects with which to 
study galaxy formation and evolution. In addition, clusters are the largest virialized ob- 
jects in the universe. By virtue of this fact we are able to, through high redshift samples, 
investigate the growth of large scale structure. A firm knowledge of the evolution of the 
cluster mass function would provide an enormous contribution to constraining cosmological 
parameters such as cxg (the normali zation of t he density perturbation spectrum) and w (the 
dark energy equation of state; e.g., I VoitI 120051 ). 



Two ingredients are required to achieve this goal. First, a statistically significant sample 
of clusters in multiple redshift bins is needed. Second, reliable mass estimates of the clusters 
in that sample must be obtained. Difficulties in reaching the first requirement include the 
vast amount of telescope time required to carry out such a search in the X-ray, and the 
propensity for false detections due to projection effects in optical surveys. The primary 
challenge in reaching the second part of this goal is again the high cost of observing time to 
achieve either X-ray or dynamical mass estimates. 



The Red-sequence Cluster Survey (RCS; ICladders fc Yeell2000l . l2005l : lYee et al. 1120071 ) 
has attempted to evade such difficulties - RCS is an optical survey which uses the color- 
magnitude relationship of cluster ellipticals to find galactic overdensities in small slices of 
redshift space. This technique has been estimated to bring false detection rates down to 
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~ 5 - 10% flGladders 1 12002 : Isiindert et al. II2007I : ICohn et al. Il2007h . The chosen filters {Rc 
and z'), optimize this finding algorithm for the redshift range 0.2 < z < 1.2 and provide 
photometric redshift information with accuracies of ~ 10%. In addition, optical richness 
information is immediately available from the survey data and, if sufficiently calibrated, this 
information could provide a highly efficient way to estimate the masses of cluster candidates. 



The first phase of the Red-Sequence Cluster Survey (RCS-l: lGladders fc Yedl2005l ). from 
which our cluster sample was drawn, covers 90 square degrees and was performed at CFHT 
and CTIO. RCS-1 has identified 6483 cluster candidates in the redshift ran ge 0.2 < z < 1.2, 
oyer 1000 of which are at least as optically rich as Abell class clusters (iGladders &: Yee 
2005h . 



The motivations for this work are to characterize high-redshift optically-selected cluster 
samples, probe cluster evolution, and move forward in attempts to calibrate a robust rela- 
tionship between optical richness and cluster mass. This paper presents a detailed analysis 
of the Chandra data of thirteen RCS-1 clusters with redshifts in the range 0.6 < 2 < 1.1 
Our analysis investigates the temperatures and gas distributions of ten of these clusters, and 
provides mass estimates for use in the calibration of relationships between optical richness 
and cluster mass. 

We also use our results to investigate the X-ray scaling laws of our sample, and thereby 
probe redshift evolution in these relationships. To facilitate comparisons between the RCS 
clusters and lower-redshift X-ray selected samples, we make use of our previous Chandra 
analysis of the Canadian Network for Observ ational Cosrnology (CNOC) subsample of the 
Extended Mediurn Sensi tivity Survey (EMSS) ( Hicks et al. II2OO6 : Yee. EUingson. fc Carlberg 
19961 : iGioia et al. Ill990l ). This sample, with redshif ts in the range 0.1 < z < 0.6, was chosen 



based on X-ray luminosity {Lx > 2 x lO^'' erg s ^; iGioia et al. Ill990l ). 



This paper is organized as follows: In Sections [2] and [3l we introduce our sample and 
describe the basic properties of our data. In Sections H] and O we investigate gas distributions 
and obtain cluster temperatures. We derive masses for our sample in Section [61 High-z 
X-ray scaling relationships are examined in Section [TJ while correlations between optical 
richness and cluster X-ray properties are explored in Section [Si In Section M we investigate 
possible sources of bias in cluster sample selection. A summary and discussion of our results 
is presented in Section [TOl Unless otherwise noted, this paper assumes a cosmology of 
Ho = 70 km s-i Mpc"\ = 0.3, and = 0.7. All errors are quoted at 68% confidence 
levels. 
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2. Cluster Sample & Observations 

Chandra Advanced CCD Imaging Spectrometer (ACIS) observations of thirteen RCS 
clusters were taken during the period 10 April 2002 - 23 October 2005. Twelve of these 
clusters were observed with the ACIS-S CCD array, and one was observed with ACIS-I, 
with an overall range in individual exposures of 10 — 100 kiloseconds. Seven of the clusters 
in this sample were observed on multiple occasions. All multiple observations were merged 
for imaging analysis to provide higher overall signal-to-noise ratios. Those clusters with 
Atobs < 3 months between observations were merged for spectral analysis as well. Each of 
the observations analyzed in this study possesses a focal plane temperature of —120° C. 

Aspect solutions were examined for irregularities and none was found. Background 
contamination due to charged particle flares was reduced by removing time intervals during 
which the background rate exceeded the average background rate by more than 20%. The 
quiescent background was further reduced by using VFAINT mode. The event files were 
then filtered on standard grades and bad pixels were removed. Table [T] provides a list of 
each of the clusters in our sample, including their precise designation (hereafter shortened 
for simplicity), redshift, obsid, and corrected exposure information for each observation. 

After the initial cleaning of each data set, 0.3-7.0 keV images, instrument maps, and 
exposure maps were created using the CIAO 3.3.0.1 tool MERGE_ALL. Data with energies 
below 0.3 keV and above 7.0 keV were excluded due to uncertainties in the ACIS calibration 
and background contamination, respectively. Ideally, all data below 0.6 keV would have 
been excluded to minimize low energy uncertainties; however, the combined faintness and 
high redshifts of our objects require the utilization of lower energy photons as well. Point 
source detection was performed by running the tools WTRANSFORM and WRECON on 
the flux images. 

Figured] contains smoothed 0.3-7.0 keV Chandrafinx images of each of the clusters in our 
sample (produced by the CIAO tool CSMOOTH), including a comb ined image of the thre e 



objects at 2 ~ 0.9 which belong to a supercluster in the 23h field (IGilbank et al. Il2007bl ). 
As seen in the figure, this sample covers a wide range of cluster X-ray morphology, from 
very regular objects (e.g., RCS1419-I-5326), through well-detected clusters with significant 
substructure (e.g., RCS2318-I-0034), all the way to very disturbed systems (e.g., RCS2112- 
6326). It is worth noting that in Figured], the brightest part of RCS2318-I-0034 does not seem 
to lie at the center of the cluster, indicating that this object may have recently undergone 
a merger, or could at least possess an appreciable amount of substructure. Together, these 
clusters represent an assembly of some of the richest high-z (0.6 < z < 1.1) clusters in the 
RCS-1 survey. 
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3. Signal-to-Noise Ratios and Cluster Positions 



To estimate the significance of RCS cluster detections in the X-ray, we made use of 
relatively simple statistics. Counts were extracted from a 500 h^Q kpc radius region around 
the aimpoint of each observation in the 0.3-7.0 keV band (C), and also from a region far 
away from the aimpoint on the same chip which served as a background {B). Obvious 
point sources were removed from each region. Signal-to-noise ratios were calculated based 
on dividing net counts, N = C — B, hj the standard deviation, a = \/C + B. Using this 
method, twelve cluster signals were detected at a signal-to-noise ratio greater than 3, with 
the remaining object detected at S/N= 1.1 (Table [2]). 

Using adaptively smoothed 0.3-7.0 keV flux images (Figured]), we determined the loca- 
tion of the X-ray emission peak of each cluster. The images of RCS2112-6326 and RCS2156- 
0448 appear to contain multiple regions of extended emission, therefore we cannot determine 
a precise X-ray position for these objects. In the case of RCS1326-I-2903, two RCS 13h clus- 
ters lie in the field of view. The original observation was designed to observe a 2; = 1.01 
cluster at an RA, Dec of 13:26:29, +29:03:06 (J2000). Our astrometry indicates, however, 
that we are most likely detecting the emission of a lower redshift RCS cluster {z = 0.75) 
at an optical position of 13:26:31, +29:03:12. Because of the uncertainty surrounding this 
detection, we have carried both possibilities throughout much of our analysis; however, we 
include the more likely candidate (at z = 0.75) in our subsequent fitting and plots. All other 
clusters (with the exceptions of RCS21 12-6326 and RCS2156-0448) were found within 31" 
of their optical positions. Table [2] lists optical positions. X-ray positions, net counts within 
500 hjQ kpc, and signal-to-noise ratios derived from the method described above. 



A radial surface brightness profile was computed over the range 0.3-7.0 keV in circular 
annuli for each cluster. These profiles were then fit with j3 models: 



where Ib is a constant representing the surface brightness contribution of the background, 
/o is the normalization and rc is the core radius. The parameters of the best fitting models 
of the ten clusters for which surface brightness fitting was possible are shown in Table [3l 
and images of these fits are given in Figure [2l Though many of the clusters exhibit hints 
of substructure, most were reasonably well fit by a model (see Table [3] for goodness of fit 



4. Surface Brightness 




(1) 
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data). Other than somewhat low normahzations, the resuhs of surface brightness fitting are 
unremarkable (0.51 < (3 < 0.72), except in the two cases of RCS1326+2903 {(3 = 1.04), which 
lies at the edge of our detection threshold (Table [2]), and RCS2318+0034 which appears to 
not be completely relaxed, as its brightest emission is slightly offset from the center of its 
extended emission (Figured]). 



5. Spectral Analysis 



5.1. Integrated Spectral Fits and R 



2500 



Spectra were extracted from each point-source-removed event file in a circular region 
with a 300 hf^^ kpc radius. In the cases of RCS0224-0002 and RCS1419+5326 (AUs > 2 
yr), individual spectra were extracted frorn each obsid and fit simultaneously. The spectra 
were analyzed with XSPEC (lArnaudl 119961 ). using weighted response matrices (RMFs) and 
effective area files (ARFs) generated with the CIAO tool SPECEXTRACT and CALDB 
3.2.2. Background spectra were extracted from the aimpoint chip as far away from the 
aimpoint as possible. 

Spectra were fitted with single temperature spectral models, inclusive of foreground 
absorp tion. Each spectrum wa s fit with the absorbing column frozen at its measured 
value (jPickev fc LockmanI 119901 ). Metal abundances were initially fixed at a value of 0.3 
solar ( Edge fc Stewart[ 1991 ). Data with energies below 0.3 keV and above 7.0 keV were 
excluded from the fits. 



Three of the clusters did not possess enough counts to constrain a spectral fit. The 
results of the ten successful fits, combined with best fitting P model parameters from Sec- 
tion m were then used to estimate the value of R2500 f or each cluste r. This is accomplished 
by combining the equation for total gravitating mass (jSarazinlll988l ) 



Mtof (< r) 



kT{r)r f d In p ^ d In T 



Gfinip \d In r 9 In r 



(2) 



(where prup is the mean mass per particle) with the definition of mass overdensity 



(3) 



where z is the cluster redshift, and A is the factor by which the density at ta exceeds 
Pc{z), the critical density at z. Here Pc{z) is given by Pc{z) = 3H{zY /8ttG = 3HqE^/8hG, 
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where Ez = [f2m(l + zY + ^a]^/^. These equations are then combined with the density 
profile imphed from the (3 model (assuming hydrostatic equilibrium, spherical symmetry, 
and isothermality) 



-3/3/2 



Pgas{r) = Po 



1 + 



(4) 



resulting in the equation 



3/3kT 
_Gfimp{A/3)TTpc 



(5) 



(Ettori 


2000; 


Ettori et al. 


2004b) 



After the initial estimation of -R2500; additional spectra were extracted from within that 
radius, and spectral fitting was performed again. This procedure was repeated until temper- 
atures and values of -R2500 were consistent for a given spectrum. Where statistically possible, 
additional fits were performed allowing abundances to vary also. Redshifts were also fit (al- 
lowing only z, Tx, and normalization to vary) for the three clusters in our sample that do not 
have spectroscopic redshifts. We were unable to constrain a redshift for RCS1326-I-2903. The 
fits of the other two clusters resuhed mz = 0.62 ±0.01 for RCS1419+5326, and z = O.rstofs 
for RCS2318+0034, within 10% and 14% (re spectively) of photome tric redshift estimates ob- 
tained using the color of the red sequence (iGladders fc Yeell2005l ). Fits with redshift fixed 
at these values were used in subsequent analysis (Table Hj). The small uncertainties in these 
values (< 10%) do not substantially affect our analysis. 

Unabsorbed 2-10 keV luminosities within R2500 were calculated using fixed abundance 
fits. These were then converted to bolometric luminosities by scaling, using a thermal emis- 
sion model in PIMMS. For the three clusters for which spectral temperature fitting was 
impossible (RCS1417+5305, RCS2112-6326, and RCS2157-0448), spectra were extracted 
within 500 h^Q^ kpc radii, and fit in XSPEC with temperatures fixed at 4 keV (slightly lower 
than the average Tx of the sample) to determine their luminosities. Temperature uncertain- 
ties of ±2 keV were folded into the errors of these estimates. To estimate Lx(A = 500), we 
extracted spectra from within that radius for the ten clusters in our detailed analysis sample 
and again used fixed abundance fits, with temperatures also fixed at the R2500 value. The 
results of spectral fitting are shown in Table HI along with 68% confidence ranges. Bolometric 
X-ray luminosities are listed with richness measurements in Table [51 
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5.2. Tx-cy Comparisons 



Velocit y dispersions for three o f the clusters in this sample were obtained from lGilbank et al 



(l2007af) andlGilbank et al. I (l2007bl ) , and are listed in Table [61 Using the a — Tx relationship 



of IXue fc Wul J2000h : a = 10^'^^ T^^^, we find that our temperatures are in agreement with 



the clusters' measured velocity dispersions in all cases (Table [H]). This result indicates that 
these three systems, at least, are not overly disturbed. 



6. Mass Estimates 

An isothermal cluster whose surface brightness is well fit by a model can be shown to 
have a gas density profile which follows Equation m Using this relationship and the equation 
of hydrostatic equilibrium (Equation [2]), total mass can be determined via 



tot 



< r) 



SP kTr {r/vcY 



G unip 1 + (r/rc 



(6) 



To estimate gas mass (again assuming hydrostatic equilibrum, isothermality and spheric- 
ity), the first step is to obtain a central density {hq = pQ/nip). There are two complementary 
ways to go about this. One is to use the surface brightness normalization: 



no 



1/2 



(7) 



where the F function results from surface brightness integration, (3 comes from the fit to 
surface brightness, He is the mean atomic mass per free electron (0.62), Xh is the hydro- 
gen mass fraction (0.707), Sq is the gas emissivity, Iq is the best fitting surface brightness 
normalization (corrected for absorption), and is the core radius. 

A second method of estimating central density makes use of both imaging and spectral 
fitting: 



"'ang 



'1 + zY K 



0.82 47rr3 EI 



cm 



Here K is the normalization of the XSPEC model and EI is the emission integral, estimated 
by integrating the (spherical) emissi on from the source out to som e radius - in our case we 
use 10 Mpc following the method of lEttori. Tozzi &: Rosati J2003h . 
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For the RCS sample we employed both of these methods, as it was crucial to confirm that 
we were not underestimating central density in these comparatively low luminosity objects. 
We also added the data from our previous Chandra analysis of the moderate redshift CNOC 
sample, to cover a wider range of redshifts in our comparison. We found that the methods 
agree (on average) to within 10%, and proceeded in our analysis using the surface brightness 
normalization method. 

From these equations, along with Equations S] and and using the results of spectral 
and surface brightness fitting, gas masses and total masses were determined out to R2500 
and R500 for the clusters in this sample. We also calculate gas mass fractions for the RCS 
clusters and find them to be systematically lower than the gas fractions of lower redshift 
X-ray selected clusters. The robustness and implications of this result are explored in detail 
in Section [731 

Gas masses, total masses, and gas mass fractions can be found in Tables [7] and [81 We note 
that while extrapolations to larger radii are possible using out measured /3-fit parameters, 
R2500 is the radius to which we have confident measures for our entire sample. 



7. Cluster Scaling Relations 



Studying the relationships between global cluster properties (Lx, Tx, Mtot, etc.) over a 
broad range in redshift allows us to investigate the infiuence of non-gravitational processes 
on cluster formation and evolution. On a less grand scale, these relationships can also lead 
to interesting clues regarding an individual cluster's dynamical state and composition, as 
well as provide a method of comparison between different cluster samples. In this paper we 
investigate the evolution of scaling relationships over the redshift range 0.1 < z < 1.0, and 
use them to characterize high-z optically-selected RCS clusters (0.6 < z < 1.0) 

To facilitate comparisons between the RCS clusters and lower-redshift X-ray selected 
samples, we make use of our previous Chandra analysis of the Canadian Network for Ob- 
servation al Cosmology (CNOC) subsample of the Extended Medium Sensi tivity Survey 
(EMSS) Jmcks et al. IboodlYee. EUingson. fc Carlberg||l996l : lGioia et al. Ill990h . The CNOC 



sample was assembled prim arily based on X-ra y luminosity, with a cut at 2 x 10 erg s in 



the original EMSS catalogs (ICioia et al. Ill990l ). and covers a redshift range of 0.1 < z < 0.6. 



This sample is not well-matched in redshift to our RCS clusters, but it is one of the 
best-studied moderate-redshift cluster samples today, with substantial information about 
both X-ray and opt ical properties ava ilable. Our previous analysis of this sample using the 
same methodology (iHicks et al. II2006I ) also allows us to make a confident comparison of our 
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measurements. 

All relationships {Lx-Tx, L^— Mtot, Mtot — Tx, Lx-Yx, and Mtot — Yx) are fit within ei- 
ther R2500 or R-500! have been scaled by the cosmological factor = H{z)/Hq = [^2^(1 + z)^ + ^a]^ 
and are fit with the form 



log 



10 Y = CI + C2 logio X. 



(9) 



In all relationships Tx is in units of 5 keV, Lx in units of 10^^ erg s ^, total mass in 10^^ M©, 
and Yxi= MgTx) in 4 x 10^^ M© keV. Best fitti ng relationships are determined using the 



bisector modification of the BCES algorithm in lAkritas fc Bershadyl (119961 ). and we cal- 



culate scatter along the Y-axis as [ S^-i^Af (login Y,- — — C2 logio /A^] , facilitating 



comparisons to previous work (e.g.. lEttori et al. Il2004al ). 



In all fits the cluster RCS0439-2904 was left out, due to spectroscopic indications (iGilbank et al 



2007al : ICain et al. 1 120071 ) that it does not consist of a single virialized mass but two closely 
spaced objects in projection along the line of sight. We perform fits at A = 2500 on the in- 
dividual samples; RCS (0.62 < z < 0.91) and CNOC (0.17 < z < 0.55), as well as combined 
data from all 23 clusters (0.17 < z < 0.91). At A = 500 only the RCS data are fit. All fits 
are then reproduced with the slope fixed at the expected self-similar value. The following 
discussions pertain to fits with two free parameters unless otherwise noted. Results from the 
fitting performed in this section can be found in Tabled while Table [TO] provides comparison 
fits from the literature. 



7.1. The Lx - Tx Relationship 



In the absence of significant preheating and/or cooling, theory predicts that cluster 
luminosities should scale as Lboi oc T^. However, observational studies have resulted in rela- 
tionships which fall closer to Lhni oc T^ (jWhite. Jones fc Formanlll997l : lAUen fc Fabianlll998 



MarkevitchI Il998l : lArnaud fc EvrardI Il999l ). These departures from theoretically expected 
self-similar scaling la ws indicate the effects of non-gravitational processes, such as galaxy 
formatio n (IVoitll2004f). There is also interest in whether the Lx-Tx relationship evolves with 
redshift (lEttori et al. Il2004al ). which we investigate in this section along with the properties 
of our sample. 

The best fitting relationships and their scatter are given in Table [9] and are plotted in 
Figure [31 At both radii (R2500 and R500), the slope of the RCS fit is found to be consistent 
with a predicted self-similar slope of 2. The CNOC sample, with a slope of 2.31 ± 0.31, 
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is only marginally consistent with predicted scaling, but does agree with other low-redshift 
Lx-Tx relationships. The main difference between the RCS and CNOC fits, however, is their 
normalization, which is significantly lower in the case of the RCS fit, translating into 2.3 ± 
0.3 X 10^^ erg s~^ at 5 keV, compared to S.SIqq x 10^^ erg s~^ for CNOC. Fits to the combined 
sample have significantly higher scatter and a much larger slope (2.90 ± 0.35), inconsistent 
with self-similar evolution. This is an interesting result - taken with the individual fits it 
suggests that redshift evolution in the normalization of Lx-Tx could be perceived as evolution 
in its slope, if the fitted sample covered a broad enough range of redshifts. 

This speculation naturally leads us to the important question of whether this trend to- 
wards lower luminosity is due to the different selection of the RCS clusters, or to a general evo- 
lutionary trend with redshift. In Table [TD] we list our fit parameters along with others taken 
from the literature. Since a number of these studies use measure ments at Rnnni we compare 



our fit s using estimates extrapolated to this radius. The results of lAUen. Schmidt, fc Fabian 



( I2OOII ). from a selected sample of 2; ~ 1 clusters, are consistent in slope with our individual 
sample fits, but even higher in normalization than the CNOC fit. This is perhaps not sur- 
prising given that the clusters in their sample are relaxed lensing clusters, many of which 
have strong cooling cores which can significantly increase the central cluster luminosity. In 
general, we do not here have enough information to excise cooling cores from the RCS data; 
however,, in the CNOC dat a, we did attempt to remove these features from the cluster tem- 



peratures and luminosities (IHicks et al. II2006I ). Discrepancies between CNOC and RCS are 



thus not likely to be due to a higher incidence of cooling cores in the lower redshift sample. 



We also compare our sample to the 0.4 < 2 < 1.3 X-ray selected sample in lEttori et al. 



( l2004al ). which consists of 28 clusters taken from the Chandra archive. At A = 500, they 
find that the slope of the Lx-Tx relationship is much steeper than that predicted by self- 
similar scaling (slope=3.72 ± 0.47; Figure [3]), suggesting a negative redshift evolution in 
the relationship (i.e., cluster s at hig h- 2: have lower Lx for a given Tx)- A similar result 



is suggested by lEttori et al. I (l2004bl ) , based on simulations guided in part by low-redshift 
observations. 

Fi gure B] shows the ext rapolation of our data to Rsoo- Seven of our nine objects lie 



on the lEttori et al. I (l2004al ) relationship, suggesting at least some agreement between the 
properties of their X-ray and our optically-selected samples. We note that in general their 
higher redshift clusters also trend towards lower luminosities. While the slope of our fit to 
the RCS sample at A = 500 is inconsistent with theirs, the slope of our combined A = 2500 
sample is in agreement with their slope, another indication that we may be resolving evolution 
in slope into changes in the normalization of the relationship with redshift. Our scatter for the 
individual fits is significantly lower than theirs (aiogv < 0.20 vs. aiogY = 0.35), whereas our 
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scatter for the combined fits becomes more comparable (0.28). Thus, the RCS high-redshift 
sample appears to be at least qualitatively similar to this high-redshift X-ray selected sample, 
in support of a trend for samples of clusters at high redshift to have have lower luminosities 
at a given temperature. 



7.2. The Lx-Mtot Relationship 



Upon examining the cosmologically corrected Lx— Mtot relationship, we again see dispar- 
ity between the normalizations of the CNOC and RCS fits. This finding provides additional 
evidence that there is less gas for a given total mass in our high redshift sample. Individual 
slopes at R2500 agree with the self-similar value of 1.33, while the slope of the combined 
sample fit is higher and inconsistent with that value (Tab le M Figure H I) . Our RCS fit at 
R500 is again consistent in normalization but not slope with lEttori et al. I (j2004al ). Likewise, 



again our combined (R2500) sample slope (1.77 ± 0.15) agrees well with theirs (1.88 ± 0.42). 
Our scatter (0.16 < criogv < 0.33) is lower in all cases. 



7.3. The Mtot - Tx Relationship 



The Mtot — relationship is by far the lowest-scatter (aiogv < 0.10) relationship in this 
work (Figure[5]), though this is largely because of the degeneracy between the two parameters, 
with much of the scatter arising from differences in the spatial distribution of gas. All of 
our fits at R2500 have consistent normalizations, and all but the CNOC fit have slopes which 
agree well with self-similar predictions. Though the CNOC slope i s higher (1.83 ±0.13), it is 
in good agreement with both of our listed R 25nn comparison fits (lAUen. Schmidt, fc Fabian 
2OOII : lArnaud. Pointecouteau. fc Prattll2005l ). as are the fits of both the RCS sample and the 
combined sample (Table [TUl) . 



AtR, 



'500 



the RCS fit is in ag r eemen t with lFinoguenov. Reiprich fc Bohringeii (120011 ). lArnaud. Pointecoui 



( I2OO5I ). and lKotov fc Vikhliniru (120051 ). Our norr nalization is somewha t higher, however, than 
all of theirs, and is in disagreement with that of lEttori et al. I (j2004al ) (Figure Ej). The three 
objects that are most responsible for driving up the normalization all have gas distribu- 
tions which appear to be more concentrated than average (/? > 0.72; Table [3]), which would 
tend to drive up the total mass at higher radius, noting again that at R500 our masses are 
extrapolations. In addition, the two most outlying points consist of our least massive clus- 
ter (RCS1326+2903, also our weakest detection in the detailed analysis sample), and our 
most massive cluster (RCS2318+0034). It is worth mentioning again that in Figured], the 
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brightest part of RCS23 18+0034 does not seem to lie at the center of the cluster's extended 
emission, indicating that this object may have recently undergone a merger, or could at least 
possess an appreciable amount of substructure. 



7.4. Gas Mass Fractions 



In Section[6l we estimate the core (R2500) gas mass fractions of our high-z sample, finding 
values which are significantly lower than both the gas fractions of lower redshift X-ray selected 
clust ers, and the expected universal gas fraction {VLb/^m = 0.175 ± 0.012; ISpergel et al. 



20071 ). Taking a weighted average over our objects results in a gas mass fraction of 4.5±0.2%, 
in comparison with the CNOC weighted mean of 9.8 =fc 0.3% and values of ~ 9% found 
within R2500 in clusters with Tx > 5 keV (IVikhlinin et al. 1120061). Poor clusters and groups 



are often found to ha ve lower gas mass fractions (jDelF Antonio. Geller. fc FabricantI Il995 



Sanderson et al. II2003I ). therefore we may expect this result for the lower temperature objects 
in our sample. This, however, does not explain our findings for the higher temperature 
objects. 

To investigate further, we first performed a K-S test on the fg values of subsets of both 
samples, choosing the 8 objects in the RCS sample and the 9 objects in the CNOC sample 
with temperatures between 3.5 and 8 keV. This test resulted in D=0.875 and P=0.002, 
indicating that the gas mass fractions of the two samples are different at a confidence level 
of greater than 99%. A histogram showing the fg distributions of these subsamples is shown 
in Figure [61 

We examined the robustnesss of this result by repeating the K-S test after attempting 
to remove the effects of any possible trend in gas fraction with temperature. To do so 
we assumed that the RCS and CNOC samples can be combined and that the resulting 
apparent trend of /gas with temper ature is physical (note that the resulting relation is much 



steeper than the one suggested by lVikhlinin et al. II2006I . ; Figure E]). Under these extreme 



assumptions, the KS-test yields D=0.764 and P=0.007, thus demonstrating the robustness 
of our earlier results. 

Low gas mass fractions have p reviously been obs erved in clusters at high red s hift b y 
both the XMM-Newton n project JSadat et al. Il2005h and Ihubin. Oke fc PostmanI J2002h 
and have been predicted i n simulations of high redshift objects dNaeai. Kraytsov fc Vikhlinin 
2OO7I : lEttori et alHbood : iKravtsov. Nagai fc VikhlininI I2OO5I: lEttori et al. Il2004bl l In addi- 
tion, an SZ/WMAP study performed by lAfshordi et al. reports that ~ 35% of ex- 
pected baryonic mass is missing from the hot ICM in their 193 clusters. Redshift evolution. 
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however, may not be the only possibihty. Multiple studies have confirmed that at least some 
fract i on of their optically selected clusters have lower than expected Ly (e.g., iBower et al. 



1994 : iDonahue et al.l l2002l : iGilbank et al.l |2004 IPopesso et al. 1120071 ). therefore sample se- 
lection may also contribute to this effect. We explore selection biases in more depth in 
Section [91 

Given the possibilities present in the literature and a current lack of sufficient data to 
perform direct comparisons with significant samples matched in both mass and redshift, it is 
difficult to determine conclusively that the low gas fractions measured here are the result of 
cluster evolution. Possible physical explanations for lower gas fractions are that our clusters 



have a comparatively higher amount of ba ryonic matter in the form of stars (jVikhlinin et al. 



20061 : iNagai. Kravtsov &: VikhlininI 120071 ). that gas is still infalling (i.e., in the process of 



virialization; IPopesso et al. 1120071 ). or that some mechanism has injected e x cess e nergy into 
the gas (i.e., galaxy formation, mergers, AGN, radio jets; iNulsen et al. I (120051 ) ). thereby 



raising its entropy at high-^. Many of these processes occur wit h relatively higher frequency 
at high-redshift (e.g.. lLacey fc Colelll993l : Eastman et al. 1120071 ). thus a general trend toward 
lower gas fractions might easily be expected in high-z clusters 



7.5. Cluster Entropy 



Cluster entropy can be used as a tool for investigating the energy budget of baryons 
in clusters ( iPonman. Cannon fc Navarro! Il999l ) . Because it may provide insight into fgas 
discrepancies, we investigate it here for our two samples. The measurable quantity S = 
Tx/nl^^ can be related to thermo dynamic entropy by K = logS*. Th e canonical radius for 



measuring this quantity is O.IR200 ( jPonman. Cannon fc Navarrdll999l ). so that is the radius 
at which we present it here. 

Figure [7] shows a plot of cosmologically corrected entropy (Et^^S) vs. temperature, with 
the relationship of iPonman. Sanderson &: Finoguenovl (120031 ) overlayed {S ~ 120 T^^^ keV 
cm^). The specific entropy of the RCS clusters seems overall to be slightly higher for a given 
temperature than that of the CNOC sample. A K-S test on the clusters with 3.5 < T < 8.0 
results in D=0.431 and P=0.208, indicating a difference between the samples at an ~ 80% 
confidence level. 

As in the case of the gas mass fractions, we perform an additional K-S test after at- 
tempti ng to remove the trend in cluster entrop y with temperature. Using the relation- 
ship of IPonman. Sanderson fc Finoguenovl (120031 ) (above), the KS-test yields D=0.764 and 
P=0.007, indicating a systematic difference in the entropies of the two samples at a > 99% 
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level. Weighted means of the corrected (5 keV) entropies of the two K-S sample subsets 
result in Scnoc = 297 ± 9 keV cm^ and Srcs = 425 ± 18 keV cm^, with the RCS clusters 
having higher entropy on average by a factor of 1.43. Because fg is proportional to gas 
density, at a constant temperature fg oc 5*^/^, therefore differences in entropy between the 
two samples can account for roughly 85% of their fg discrepancy, indicating that additional 
factors may be in effect as well. 

It remains difficult to determine the relative contributions of evolution and selection to 
po ssible differences in entropy. Expectations of higher merger and AGN activity at high- 



er (ILacey &: Cold Il993l : Eastman et al. 1 120071 ) suggest that an evolutionary explanation is 



feasible; however. X-ray surveys that select high central density objects may be prone to 
preferentially pick out low-entropy systems. 



7.6. Yx Relationships 



The product of cluster temperature and gas mass, Yx = MgTx has been shown to be a 
reli able, low-scatter proxy for total clus t er mass and to be well correlated to X-ray luminos- 
ity (IKravtsov. Vikhlinin &: Nagail l2006l : iMaugharu 120071 ) . Here we investigate relationships 
between these quantities a nd Yx for our hi gh r edshift optically selected sample. W e adopt 
self- similar scaling from lMaughanl (120071 ) and lKravtsov. Vikhlinin fc Nagail (120061 ) for the 
Lx and Mass relationships, respectively; and use their best fitting slopes for our constrained 
slope fits. 

Our individual samples can be seen to lie again on two separate relationships in the 
Yx-Lx plane, of similar slope and differing normalization (Figure [S]). This can once more be 
explained as stemming from systematically lower gas naass fr actions in the RCS sample. We 
cannot make normalization comparisons to iMaughanl (120071 ) at R2500; however, none of the 
slopes of our Lx-Yx fits are consistent with the slope resulting from fits to his overall sample, 
and when we fix the slope to his value our scatter increases by a factor of ~ 2. At R500 our 
fit to the RCS data agrees neither in slope nor in normalization with his fit (Figure [S]), but 
again it should be mentioned that our R500 values have been extrapolated from data within 

R2500- 

There is an even more significant discrepancy between the normalizations of the CNOC 
and RCS samples in the Yx—Mtot relationship. This is easily explained as we are already 
aware that gas mass fractions are lower in the RCS sample, and total mass vs. Yx (oc Mgas) 
highlights this differenc e. We find overall closer agre e ment with the slope of the Yx— Mtot 
relationship modeled by iKravtsov. Vikhlinin fc Nagail (120061 ) than we did in the case of Lx- 
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Yx- Here we see marginal agreement at R2500 between their R500 slope and that of the RCS 
fit, and consistency with the slope of the CNOC sample relationship. At R500, though the 
slope of the RCS fit is still consistent with theirs, the normalizations diagree. The reason for 
this is illustrated nicely in the right panel of Figure M The four clusters in our sample which 
do not lie on their relationship are those with the lowest gas mass fractions. And again we 
see that the three biggest outliers are those with the highest /3 values, and that of these the 
two most discrepant are RCS1326+2903 and RCS2318+0034. 



8. Correlations with Optical Richness 



Optical richness is effectively a measurement of galaxy overdensity within a given aper- 
ture, normalized for the evolving galaxy luminosity function and the exp ected spatial distri- 
butio n of galaxies in the cluster. Our chosen richness measuremen t. Bg^^ (lYee fc Lopez-Cruz 
1999 ). represents the galaxy- cluster spatial covariance amplitude ( Longair fc Seldnei 19791 ). 



To 



B r"^ 



(10) 



In practice, Bgc is based on the excess number counts of galaxies within 357 kpc of 
the cluster optical center, with a normalization applied to correct for the expected spatial 
distribution (here we assume 7=1.8, which is in general agreement with actual galaxy distri- 
butions at these radii) and for the evolving luminosity function of cluster galaxies. Though 
some uncertainties exist in the evolution of cluster galaxies at redshifts oi z > 0.5, they 
can be minimized by employing a red-sequence optical richness, Bgc,red, which is calculated 
using only the more uniformly evolving red galaxies in a cluster, and which may be better 
correlated with the underlying cluster mass. This is the optical richness parameter which 
will be used throughout this work. Values of Bgc,red for this sample are given in Table [51 

I t has been showri that Bgc cor relates well w i th the X-ray parameters of relaxed clus- 
ters (lYee &: EUingsoru |2003| ). and in iHicks et al. I J2006h we have derived relationships for 
correlations of X-ray properties with Bgc^red- These relationships, however, were calibrated 
for X-ray selected clusters at moderate redshift, and therefore may not accurately describe 
our current sample. Here we test these correlations for optically selected clusters at high red- 
shift. In the following we will assume that Bg^^red behaves similar to the X-ray temperature 
when comparing to the X-ray properties. The rationale for this choice is that for a cluster 
with a fixed density profile, both the temperature and Bgc,red do not change with redshift, 
whereas for instance M2500, Lx^ etc. do change (following self-similar evolution). 
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Our actual data do not extend much beyond -R2500; so most of our fitting is performed 
within that radius. Cluster proper ties included in our fi t s are Ly, Tx, and total mass. Again 
we employ the BCES algorithm of lAkritas &: Bershadyl (119961 ). For each of our fits we adopt 
the form 



logio Y = Ci + C2 logio Bgcred 



where Y represents the particular property being fit. For Lx, Tx, and total mass, units 
of 10"^^ erg s~^, 5 keV, and 10^"^ M© were used, respectively. RCS0439-2904 was again 



removed from fitting procedures, as i t has been confirmed to be two closely-sp aced systems 
in projection along the hne of sight (IGilbank et al. I l2007al : ICain et al. 1120071 ). Best fitting 



parameters and scatters are given in Table [TTl along with comparison fits from the literature. 

Figure [To] shows the relationship between temperature and richness. Here, there is little 
evidence for a systematic difference between the samples, with both the CNOC and the com- 



bined sample showing statistical agreement with the expected slope of 2/7=1.11 (lYee fc EUingson 



2OO3I ). The RCS sample is on average slightly cooler at a given Bgc,red, a tendency that might 
stem from sample selection (see Section [HI below). Note that both X-ray and optically- 
selected samples contain a few outliers, scattering towards higher temperature or lower rich- 
ness. Figure [11] shows relationships between richness, M2500 and M2oo- At R2500 only the 
CNOC sample shows agreement with the expected sl ope of 8/7. AH R?.nn fits are consis- 



tent with the values obtained for t he CNOC sampl e bylYee fc EUingson] (120031 ) using galaxy 



dynamics, and the fit reported in ]Blindert et al. I (12007] ) for a sample of 33 RCS clusters 



(0.2 < z < 0.5). The consistency between these fits indicates a general agreement between 
both the samples and the different mass estimators; however, it should also be noted that 
error bars on the fit parameters are quite large for this relationship. 

The Lx-Bgc,red plot (Figure ]T2|) in contrast, shows quite a bit of scatter for both samples, 
and a significant offset between the RCS and CNOC samples. This offset is expected for 
the RCS sample, based on the results of Section [TTTl but here we also include the additional 
clusters with low X-ray luminosity for which Tx could not be derived. RCS0439-2904, the 
object which is spectroscopically confirmed to be a projection of two less massive systems, 
is the cross second from the right. The significantly higher amount of scatter that we see in 
this relationship when compared to any of the other Lx relationships suggests that Bgc,red is 
a less reliable predictor of X-ray luminosity than Tx, total mass, or Yx- 

Of all the richness relationships we investigate, Bg^^red is best correlated to X-ray tem- 
perature, with an average scatter (all fits) of only aiogy ~ 0.16 for the objects with measured 
Tx (minus RCS0439-2904, which is not included in fitting). Interestingly, this is less than 
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the average scatter of our Lx-Tx relationship for the same objects (all fits; aiogv ~ 0.22), 
probably again due to the issue of missing ICM baryons in the high-z RCS sample. Since 
Tx is closely related to total mass, whereas Lx is intimately tied to the gas density, differing 
gas mass fractions will shift the normalization in Lx-Tx significantly, whereas Bgc^red may 
be fairly independent of the amount of X-ray emitting gas. 

When fitting mass to richness, overall scatters average (all fits) to (Xiogv = 0.28 for 
M2500, and 0.32 for M2oo- The scatter in the RCS sample at R200 is particularly large, due 
in part to the two objects with high (3 values (RCS1326+2903 and RCS2318+0034), and 
RCS2320-I-0033 which has a lower than expected Bgc,red- In its role as a mass estimator, 
Bgcred produces on average 0.07 — 0.19 more scatter in aiogY than the Tx-based mass proxies 
investigated here, and may suffer from a fraction of objects whose richnesses are affected by 
projection. However, the comparative speed and ease with which it can be obtained still 
recommend it as a potentially useful tool for mass estimations of large high redshift cluster 
samples. 



9. Sample Selection and Biases 



The RCS sample is effectively selected by richness ( Bgp rsn), whereas the CNOC sampl e 
was compiled from objects with high X-ray luminosity (lYee. EUingson. &: Carlberd Il996l ). 
Discrepancies in the relations between X-ray and optical properties for these two samples 
may thus partially be caused by sample selection, esp ecially if the uriderly ing distribution of 
X-ray-to-optical properties is intrinsically broad (e.g.. iGilbank et al.ll2004l ). Both X-ray and 
optical surveys would then be expected to deliver biased samples of clusters, with the degree 
of bias based on the level of scatter in the selection criterion. Here we discuss three sources 
of selection bias in optical and X-ray cluster samples: optical projection effects, Eddington 
bias due to observational uncertainty, and sample bias for both optical and X-ray samples. 

One important difference in the cluster samples stems from the RCS cluster-finding 
process. While all but one of our objects were confirmed (to S/N> 3) as extended X- 
ray sources, the RCS sample is expected to also include a small fraction of objects whose 
richness is boosted by the superposition of other structures having galaxy colors similar 
to the cluster's red sequence. While these projections are much less problematic than in 
monochromatic cluster searches, they may still add systems into the R-CS cat alog with true 
richnesses significantly lower than the measured Bgcred- iGladders I ((2002) performed a 



series of simulations which suggest that the fraction of RCS clusters composed of significant 
projections is on the order of 5-10%. This estimate has been c onfirmed at 2; ~ 0.3 via 
extensive spectroscopy of 33 RCS clusters (IBlindert et al. 1120071 ): and at z ~ 0.8 from a 
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sample of 12 clusters (iGilbank et al. Il2007al ). Additional s pectroscopy a s well as weak lensing 



estimates of additional clusters is underway. Recently, ICohn et al. I (120071 ) examined the 



effects of local structures on the red sequence using the Millennium cosmological simulations. 
They found that in the simulation, the frequency of significant projection increases at higher 
redshift, to ~ 20% at 2; = 1. However, their cluster- finding algorithm and richness estimate 
were significantly different from the RCS algorithm in many aspects (galaxy magnitudes and 
colors, radial extent and background corrections), so this may not be directly comparable to 
the samples discussed here. 

Our Chandra observations suggest that perhaps 3 of 13 observed clusters may have X-ray 
luminosities which are significantly lower than expected from the RCS Lx-Bgc,red relationship. 
All three of the outliers in Figure [12] have been observed in detail spectroscopically, and two 
of these were found to have at least some degree of overlap with additional structures in 
the hne of sight. RCS0439-2904 was found spect roscopically to consist of two objects i n 
such close proximity that they may be interacting (IGilbank et al. Il2007al : ICain et al. 1120071 ). 
RCS 141 7-1-5305 is a similar case, though here the overlapping sy stems are different enough 
in redshift that they might be unrelat ed (jcilbank et al ~ 2007a ). RCS2112-6326 exhibits a 
single spectroscopic peak at z ~ 1.1 (IBarrientos et al. 1120071 ). It is not clear whether the 
highest richness systems in this sample might be subject to a higher contamination rate than 
the RCS-1 survey as a whole. 

Because there is significant observational uncertainty in our richness estimates, and the 
number of clusters declines rapidly with increasing richness, it is also necessary to evaluate 
effects of a possible Eddington bias in the X-ray/optical relationships. We calculate this 
possible bias by using the observed distribution of Bgc,red in the RCS-1 sample, which falls 
as Bgc^red^, where N ~ —4. Uncertainties in Bgc^red are calculated based on the statistics of 
galax y counts in the clusters and in the statistical foreground/background galaxy distribu- 
tion (lYee fc Lopez-Crug|l999l ). and tend to increase modestly with increasing richness. We 
model the typical gaussian 1-a uncertainty in Bgc,red as a function of Bgc,rod from an empirical 
fit to the observed distribution in RCS-1: 



logio (^) = 0.899 + 0.535 log^o (Bg^^ed) (12) 

This relationship predicts that the uncertainty will be ~ 180 at Bgc,rcd=300 h^Q Mpc^-^^, 
at the lower end of our cluster richness distribution, and ~ 320 at 1000 h^^ Mpc^'^^ for very 
rich clusters. (Note that these errorbars are not identical to the detection significance for 
the cluster, but instead reflect the uncertainty in the measurement of the cluster's richness). 
Convolving these relations predicts that the true distribution in richness for a measured 
Bgcred is skewed to lower values, with a mean value that is ~ 80 — 90% of the measured 



- 21 - 



value. We then use our observed relationship between X-ray temperature derived from the 
CNOC clusters and B^c to calculate that the mean observed temperature for RCS clusters 
should be about 10% lower than the expected relationship at Bgc,rcd=1000 h'^Q Mpc^ ''^ and 
about 20% lower at Bgc,red=500. These decrements will also tend to steepen the logarithmic 
slope of the Tx-B^c relationship by about 0.15. Varying the distribution parameters within 
reasonable limits produces corrections on the order of 10-30% in normalization at a given 
Bgcred, and a systematic increase of 0.1-0.3 in the slope. Comparison with the T^-Bgc^red re- 
lation shown in Figure [TU] suggests that a correction for this Eddington bias would ameliorate 
the discrepancies between the RCS and CNOC samples, likely resulting in statistical agree- 
ment between their respective fits. 



A similar calculation for X-ray luminosities was performed, with decrements in the X- 
ray luminosity of about 40% and 25% beneath predicted values at Bg^^red 500 and 1000 h^^ 
Mpc^'^^, steepening the logarithmic slope of the -Lx-Bgc,rcd relationship by about 0.3. This 
correction is not, however, sufficient to create agreement between the RCS and CNOC sam- 
ples once bias in richness measurements has been accounted for, as is also indicated by their 
differing L^-Tx relations. Note these calculations assume that observational uncertainty is 
the primary source of scatter in the correlations. 

A final consideration in comparing X-ray and optically-selected samples is the possibility 
that both selection methods produce biases when selecting clusters from a population with a 
significant intrinsic variation in X-ray or optical properties. If there is a significant intrinsic 
scatter in the properties of gas in cluster cores, systematic differences in X-ray characteris- 
tics between optically and X-ray selecte d samples may natur ally arise. The ROXS survey. 



a joint X-ray/optical survey for clusters (IDonahue et al.ll2002l ) found that some of their op 



tically selected clusters had lower th an expected Lx, suggesting that selection effects could 



be culpable. iGilbank et al.l (120041 ) also performed an independent X-ray/optical survey 
for clusters, using the red-sequence as well as the monochromatic matched-filter technique. 
They found that the red-sequence methodology significantly out-performs monochromatic 
techniques in discovering and characterizing clusters. Even so, they also found a significant 
difference in the X-ray luminosities of X-ray versus optically-selected clusters, with s evera l 



examples of spectroscopically-confirmed low Lx clusters. More recently IStanek et al. I (120071 ) 
report that Malmquist bias may be responsible for a higher (by a factor of ~ 2) average Lx 
in X-ray fiux-limited samples. These studies all suggest that for a given cluster mass or 
temperature there is a significant intrinsic scatter in X-ray luminosity or optical richness, or 
possibly both. 

If this is the case, then both X-ray and optically-selected clusters may be prone to bias. 
We first consider the effects of such bias on our observed Lx-Tx relations, where we found 
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the RCS clusters to be systematically lower in luminosity for a given temperature, and in 
addition calculate lower core gas fractions. This discrepancy could be interpreted as evidence 
for evolution in the properties of the ICM, and the loose agreem ent of the RCS data with 
the high redshift X-ray selected sample of lEttori et al. I (l2004al ) supports this conclusion. 
In addition, our X-ray luminous CNOC comparison sample m ay include a signifi cant bias. 
These clusters were chosen from the wide-area EMSS survey ( iGioia et al. Ill990l ) primarily 
based on their X-ray luminosities and may indeed represent a sample of particularly luminous 
clusters. We summarize by noting that these selection biases in both X-ray and optical 
samples can be significant, but can be evaluated quantitatively, given additional independent 
information about the underlying cluster mass. In general, variations in the X-ray properties 
of clusters can be inferred most robustly from optically-selected clusters, and vice-versa, to 
minimize these effects. 



10. Summary and Discussion 

We have performed an in-depth X-ray investigation of 13 high redshift (0.6 < z < 1.1) 
optically selected clusters of galaxies from the Red-sequence Cluster Survey (RCS; Tabled]). 
All but one of these clusters was detected by Chandra at a signal-to-noise ratio of greater than 
3 (Table [2]), though two additional clusters in our sample (RCS1417+5305 and RCS2112- 
6326), did not posess enough signal to support further analysis. Initial imaging of the objects 
reveals that the RCS sample spans a wide range in cluster morphology (Figure [1]), from very 
regular objects (e.g., RCS1419+5326) to more disturbed systems (e.g., RCS2112-6326). 

Surface brightness profiles were extracted for ten clusters in 1 — 2" annular bins, and 
were reasonably well fit by single j3 models (Figure [2] and Table [2]). Cluster emission was 
modeled with XSPEC, beginning with a spectral extraction region of 300 kpc radius. The 
results of single temperature spectral fits, combined with best fit (3 models, were used to 
determine R25oo- Spectra were re-extracted from regions of that radius for further tempera- 
ture fitting and R2500 luminosity estimates, until extraction regions and R2500 estimates were 
in agreement. Results of this process are given in Table HI We have also used the a — Tx 
relationship to comp are the X-ray temperatures of three of our obj ects to currently available 



velocity dispersions (IGilbank et al. Il2007al : iGilbank et al. Il2007bl ). We find consistency in 



all cases (Table E]), suggesting that these three objects are at least relatively undisturbed. 

Using the results of both spectral fitting and surface brightness modeling. X-ray masses 
were calculated for ten clusters in our sample out to R2500 (Table [7j), with extrapolation to 
R500 (Table [8]). Canonical X-ray scaling laws were investigated for nine clusters and com- 
pared to those of the moderate redshift (0.1 < 2; < 6) CNOC sample (Tables P and [TU]). 
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For the Lx — Tx relationship (Figure [3]), both RCS and CNOC fits have slopes consistent 
with self-similar predictions; however, their normalizations disagree. Interestingly, the slope 
of our c ombined RCS-CNOC sample agrees with that of X-ray selected clusters at similar 
redshift (lEttori et al. Il2004al ). suggesting that evolution in the normalization of the Lx — Tx 
relationship may lie behind the observed steeper slopes. Results from Lx— Mtot fits are qual- 
itatively very similar to those of the Lx-Tx relationship (Figure H]), RCS and CNOC slopes 
are consistent with self-similarity at R25oo; but disagree in normalization due to differing 
ICM densities. 

The most notable outcome of our mass estimations is that the A = 2500 gas mass 
fract ions of RCS clusters are lower than expected by a factor of ~ 2 (jVikhlinin et al. 



20061 ). Low gas mass fractions are also r eported in the findings of other high redshi 



2007 



ter studies, both observation and theory (Lubin. Oke &: Postmanll2002l: lEttori et al. 



't clus- 
2004^ ^. 



Kravtsov. Naeai fc Vikhlininl boosi: ISadat et al.lboosi : lEttori et all2006l : iNaeai. Kravtsov fc Vikhhnin 



Afshordi et al. 1120071 ). Physical explanations for low gas fractions include suggestions 



that much of the baryonic mass has been converted into stars, that the gaeous component 
of these clusters is still infalling, or that some mechanism (i.e., galaxy formation, AGN, 
mergers, radio jets) is responsible for raising the entropy of the gas. Though we do see some 
evidence for higher entropy in the RCS sample (Figure [7]), it may not be enough to explain 
the entire fgas discrepancy. Further study will have to be undertaken to determine whether 
the overall lower gas fractions that permeate this sample are ubiquitous at high redshift, or 
an outcome of sample selection. 

Explanations aside, the growing evidence that massive {Tx ~ 6 keV) clusters may have 
an evolving or broad range of central gas mass fractions may have important consequences for 
the interpretation of future cluster surveys in the microwave and X-ray bands, which select 
in part on the basis of central gas density. Scatter in this parameter will tend to reduce 
completeness and, if not properly accounted for, inject bias in the cluster samples such 
surveys produce. In addition, si nce SZ mass determin ations tend to rely on the assumption 
of a constant gas mass fraction (IReid &: Spergelll2006l ). complimentary data may be required 
to avoid systematic errors in SZ total mass estimates. 

Using red-sequence optical richness measurements of both samples the relationships 
between Bgc,red and global cluster properties {Tx, Lx, M2500, and M200) were investigated 
(Table [TTI). We find that Bgc,red is poorly correlated to X-ray luminosity, with average scat- 
ter in the relationships reaching cxiogY ~ 0.36. Temperature, however, is nicely predicted 
by optical richness. The scatter in the T^— Bgc^red relationship is lower even than the scat- 
ter we find in Lx — Tx (criogY ~ 0.16 compared to 0.22), though this measurement may 
exclude several of the strongest outliers in our sample. This can be explained by the lack 
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of hot baryons in RCS clusters, which would affect cluster Lx, but not temperature or 
Bgcred- Richn ess- mass relationships are generally consiste nt with one another and with pre- 



vious studies (lYee fc Ellingsorul2003l : iBlindert et al. 1120071 ). and we find an average scatter of 



c"iogY ~ 0.30 for these relationships (Figure [TTl). Though this scatter is on average somewhat 
higher than the other mass proxies investigated here (by 0.07 — 0.19 in criogv), the compara- 
tive speed and ease with which it can be obtained still recommend it as a promising tool for 
mass estimations of large high redshift cluster samples. 
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Fig. 1. — Smoothed Flux Images. Adaptively smoothed X-ray flux images of our sample 
in the 0.3-7.0 keV band. Circles denote calculated values of R2500 for each cluster. The three 
single cluster images which lack circles did not contain enough cluster signal to constrain a /3 
model or a temperature, and thus lack estimates of R25oo- In each image, north is up and east 
is to the left. The last image shows the three clusters which make up the z = 0.9 supercluster 
in the 23h field. The aimpoint cluster (RCS2319+0038) lies at the top of the image on 
the backside illuminated CCD ACIS-S3, and the other two clusters (RCS2319+0030 and 
RCS2320+0033) lie on the frontside illuminated CCD, ACIS-S2. Instrumental differences 
in the two chips cause their respective backgrounds to have slightly different values in the 
image. 
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Fig. 2. — Surface Brightness Profiles. Radial surface brighness profiles for the 0.3-7.0 keV 

band accumulated in annular bins for ten clusters in our sample. A solid line traces the best 
fitting P model of each cluster. Horizontal dotted lines represent best fit background values, 
and vertical dashed lines indicate R25oo- Many of the profiles exhibit some substructure; 
however, most were reasonably well fit by a standard f3 model (see Table [3] for goodness of 
fit data). 
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Lx vs. Tx (A = 2500) 
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Fig. 3. — Lx-Tx Relationships. Left panel: X-ray temperatures are plotted against cos- 
mologically corrected unabsorbed bolometric luminosities within R25oo- Squares designate 
moderate redshift CNOC clusters (2;avg = 0.32), and diamonds represent high- 2; RCS clusters 
(zavg = 0.80). The dashed line traces the best fitting relationship for only the RCS clusters, 
which has a slope of 2.05 ± 0.3 and the dot-dash line denotes the best fit to the CNOC data 
with a powerlaw slope of 2.3 ± 0.3, both in agreement with self-similar expectations. The 
solid line indicates the best fitting relationship for the entire sample, wi th a slope of 2 . 9=1=0.3 , 
inconsistent with the self-similar value, but in marginal agreement with lEttori et al. I (l2004al ) 
who find 3.7 ± 0.5 for a cluster ensemble with 0.4 < z < 1.3. Right panel: Lx vs. Tx at 
A = 500. The solid line denotes our best-fitting relationship for the RCS clusters with 
sl ope 1.8 =1= 0.4, again c onsistent with self-similar scaling. The dot-dash line shows the fit 



of lEttori et al. I (j2004al ). which was also measured within Rsoo- Seven of our nine objects 
lie on their relationship, suggesting at least some agreement between the properties of their 
X-ray and our optically-selected samples. 
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Fig. 4. — Lx— Mtot Relationships. Left panel: X-ray mass is plotted against cosmologically 
corrected unabsorbed bolometric luminosity within R25oo- Diamonds designate high redshift 
RCS clusters, and squares represent moderate redshift CNOC clusters. The dashed line 
traces the best fitting relationship for only the RCS clusters, which has a slope of 1.4 ± 0.1, 
while the dot-dash line denotes the best fit to the CNOC data with a powerlaw slope of 
1.3 ± 0.2, both again in excellent agreement with the self-similar slope of 1.33. The solid 
line indicates the best fitting relationship for the combined sample, which again has a higher 
slope of 1.77 ± 0.15. Right panel: Our Lx-M^rm data i s plott ed with both our relationship 
(solid line; slope 1.03 ± 0.28) and that of Isttori et all (|2004al ) (dot-dash; slope 1.88 ± 0.42) 
overlayed. 
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OSS vs. Tx fA = 2500 
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Fig. 5. — Mtot — Tx Relationship. X-ray temperatures are plotted against cosmologically 
corrected mass estimates from Section [61 Squares designate the CNOC clusters (0.1 < z < 
0.6), and diamonds represent RCS clusters (0.6 < z < 1.0). The dashed line indicated the 
RCS fit, with a slope of 1.5±0.3. Fits to the entire sample (solid line) are also conisistent with 
self-similar evolution (slope 1.6 ± 0.2). The CNOC fit has a high er slope (dot- dash; 1.83 =fc 
0.13) , but is in good agreement with those from the literature (j Allen. Schmidt, fc Fabian 
200ll : lArnaud. Pointecouteau. fc Prattll2005l ) 
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Fig. 6. — Gas Mass Fractions. Left panel: Tx is plotted against gas mass fractions within 
^2500- Squares designate moderate redshift CNOC clusters (0.1 < z < 0.6), and diamonds 
represent higher-^ RCS clusters (0.6 < z < 1.0). The solid line indicates the best fitting 
relati onship for the entire sa mple, with a slope of 1.0 ± 0.2. Circles indicate points taken 



from Vikhlinin et al. 



J2006h . Right panel: Histogram of gas mass fractions for the eight 
RCS (left side) and nine CNOC clusters (right side) with 3.5 < Tx < 8 keV. A K-S test 
performed on these two samples resulted in D=0.875 and P=0.002, indicating that the gas 
mass fractions of the samples are different at > 99% confidence. 
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Fig. 7. — Cluster Entropy. Entropy measured at O.IR200 is plotted against X-ray temper- 
ature for the clusters in this study. RCS clusters (diamonds) have slightly higher entropies 
for a given Tx than CNOC clusters (squares), though probably not enough to account for 
the whole of the discrepanc ies seen in gas fractions between th e samples. The dashed line 
indicates the relationship of iPonman. Sanderson fc Finoguenovl (120031 ). 
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Fig. 8. — Lx-Yx Relationship. Yx is plotted against X-ray luminosity within R25oo- 
Though CNOC (squares; dash-dot line) and RCS (diamonds; dashed line) slopes agree, their 
normalizations are inconsistent, probably due to differences in gas mass fractions between 
the two samples. The solid line indicates the best fitting relationship for the entire sample. 
Overall this is the tightest relationship involving Lx that we investigate in this work. 
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Fig. 9. — Mtot — Yx Relationship. Left panel: Yx is plotted against mass estimates 
within R2500- Here we find marginal agreement between the slopes of our individual 
fit s: CNOC (squares; dot-dash line ) and RCS (diamonds; dashed line), and the slope 
of lKravtsov. Vikhlinin fc Nagail (120061 ). Right panel: The Fy— M^,nt relationship is fit for RCS 
cluste rs at A = 500 (solid line). In a direct comparison with iKravtsov. Vikhlinin fc Nagai 
( I2OO6I ) (dot-dash), our slopes are in agreement but normalizations are inconsistent. The 
points that do not lie on their relationship are also the clusters that have the lowest gas 
mass fractions in our sample. Bear in mind that we are extrapolating to get out to R5005 as 
our data mostly lie within R25oo- 
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Fig. 10. — Tx vs. Bgcred- A log-log plot of Tx vs. Bgcred witliin R2500 for the CNOC 
(squares) and RCS (diamonds) samples. Error bars represent 68% confidence intervals. The 
fit to the combined sample (solid line) is in agreement with both the CNOC (dot-dash) and 
RCS (dashed) fits. Bot h the CNOC and comb ined fits are also consistent in slope with the 
expected value of 1.11 (lYee fc EUingsonl |2003| ). These relationships, on average, show the 
lowest scatter of any richness relationships investigated in this work, with a lower average 
scatter even than Lx-Tx- The scatter of the RCS fit is mostly driven by the outlying point 
(RCS2320+0033) which has a very low Bgc,redfor its mass. 
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Fig. 11. — X-ray Mass vs. Bgc^rcd- Left panel: Total mass is plotted against Bg^^red 
for R<R25oo- Error bars represent 68% confidence intervals. CNOC clusters are shown as 
squares, and diamonds designate the RCS sample. The dashed line shows the best relation- 
ship for the RCS sample, while dot-dash and solid lines indicate fits to the CNOC and com- 
bined samples, respectively. The average scatter in the relationships is aiogv = 28%; however, 
all fits are consistent with one another due to large errors. Right panel: Bp-r. red vs. M2nn- X- 



ray ni asses were extrapolated to R200 for comparison with the relationship of lBlindert et al. 



(120071 ). All three of our fits are consistent with their relationship, which was determined via 
dynamical investigations of 33 moderate redshift RCS clusters. Scatter in our fits averages 
to 32% 
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Fig. 12. — Lx vs. Bgcrcd- A log-log plot of Lx vs. Bgc,red for the combined CNOC/RCS 
sample. CNOC clusters are shown as squares, and diamonds designate RCS clusters that are 
included in fitting. Asterixs indicate the four clusters that were not included in fitting due to 
non-detection (point farthest to the left), known superposition (RCS0439-2904; point second 
from the right), or insufficient counts (remaining two points). Luminosities for these objects 
(excluding RCS0439-2904) were determined in XSPEC assuming a temperature of 4±2 keV. 
Error bars represent 68% confidence intervals. A dashed line shows the RCS sample fit, while 
dash-dot and solid lines indicate CNOC and combined sample fits, respectively. On average, 
these relationships show the most scatter of all those that we investigate in this work. 
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Table 1. Cluster Sample 



Cluster 


z 


1" 


obsid 


Individual Exposure 
[seconds] 


Total Exposure 
[seconds] 


RCS022434-0002.5 


0.778^ 


7.44 


3181 


12051 


100844 








4987 


88793 




RCS043938-2904.7 


0.960'^ 


7.93 


3577 


64507 


93263 








4438 


28756 




RCS110723-0523.3 


0.735^ 


7.28 


5825 


49466 


94058 








5887 


44592 




RCS132631+2903.1 


0.75'^ 


7.34 


3291 


30907 


65499 








4362 


34592 




RCS141658+5305.2 


0.968= 


7.95 


3239= 


62824 


62824 


RCS141910+5326.2 


0.62f 


6.79 


3240 


9904 


57307 








5886 


47403 




RCS162009+2929.4 


0.870= 


7.71 


3241 


35953 


35953 


RCS211223-6326.0 


1.099S 


8.17 


5885 


70520 


70520 


RCS215641-0448.1 


1.080^5''' 


8.14 


5353 


36558 


71259 








5359 


34701 




RCS231831+0034.2 


0.78^ 


7.44 


4938 


50454 


50454 


RCS231953+0038.0 


0.900' 


7.79 


5750 


20902 


74539 


(RCS231948+0030.1) 


(0.904)' 


(7.80) 


7172 


17947 




(RCS232002+0033.4) 


(0.901)' 


(7.79) 


7173 


20899 










7174 


14791 





''Hicks et al. 


(2007) 




'^Cain et al. 


(2007) 



Gilbank et al. I (|2007aD 



'^From photometric data (iGladders &: Yeell2005l ). This cluster may be at 2; ~ 1.01 (see text 
for explanation). 



''ACIS-I observation. 
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^From X-ray spectra (this work), see text. 



TOarrientos et al. I (120071 ) 

^ID uncertain (see iBarrientos et al. 1120071 ). 



Gilbank et al. I (l2007bh 
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Table 2. Cluster Positions and Detection Details 



Cluster 


Optical Position^ 


X-ray 


Position^ 


Separation 


Net Counts^ 


S/N Ratio 




RA 


Dec 


RA 




Dec 


["] 






RCS0224-0002 


02:24:34.1 


-00:02:30.9 


02:24:34 


2 


-00:02:26.4 


4.7 


1102 


16.2 


RCS0439-2904 


04:39:38.0 


-29:04:55.2 


04:39:37 


6 


-29:04:50.3 


7.2 


461 


6.5 


RCS1107-0523 


11:07:23.4 


-05:23:13.7 


11:07:24 





-05:23:20.7 


11.4 


1056 


15.5 


RCS1326+2903^ 


13:26:31 


+29:03:12 


13:26:31 


3 


+29:03:31.0 


19.9 


181 


3.1 


RCS1417+5305 


14:16:59.8 


+53:05:12.2 


14:17:01 


5 


+53:05:16.2 


15.8 


138 


4.7 


RCS1419+5326 


14:19:12.1 


+53:26:11.0 


14:19:12 


1 


+53:26:11.6 


0.6 


2903 


40.2 


RCS1620+2929 


16:20:10.0 


+29:29:21.5 


16:20:10 


1 


+29:29:20.8 


1.5 


257 


7.1 


RCS2112-6326 


21:12:23.1 


-63:25:59.5 










232 


4.8 


RCS2156-0448 


21:56:41.2 


-04:48:13.3 










54 


1.1 


RCS2318+0034 


23:18:31.5 


+00:34:18.0 


23:18:30 


8 


+00:34:02.5 


19.9 


1161 


21.5 


RCS2319+0030 


23:19:48.7 


+00:30:08.5 


23:19:46 


8 


+00:30:14.3 


29.1 


780 


17.8 


RCS2319+0038 


23:19:53.9 


+00:38:11.6 


23:19:53 


2 


+00:38:12.5 


10.5 


1742 


26.2 


RCS2320+0033 


23:20:03.0 


+00:33:25.1 


23:20:02 


1 


+00:32:57.6 


30.6 


725 


16.8 



^All positions are given for equinox J2000. 
^^0.3-7.0 kcV band, within R < 500 /ifo^ kpc 

'^Tliere is also an RCS cluster at 13:26:29, +29:03:06, which is 39.2" from the X-ray centroid (see text for 
details) 
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Table 3. /3-Model Fits 



Cluster 


Tc [hjo kpc] 


/3 


-'0 


T a 
-'B 


XVDOF 


RCS0224-0002 




n 79+0.04 
' ^-0.04 


7 1 +0.4 
' •-'--0.4 


9 qq+0.02 
^•'-''-'-0.02 


210.7/199 


RCS0439-2904 


1081? 


0.591H1 


c, 9+0.4 
^■^-0.3 


q on+0-03 


67.8/64 


RCS1107-0523 


311^ 


0.511^1 


6211 


9 c:+0.1 
^•'-'-0.1 


216.1/191 


RCS1326+2903 


148+^^ 


1.041°:°^ 


o 4+0.3 
"-••^-0.2 


2 07+0-03 
-0.03 


54.9/64 


RCS1419+5326 


52+1 


0.601^} 


1891? 


9 0/1+0.02 
^■'J^-0.02 


170.6/155 


RCS1620+2929 




0.601H! 


13.3li 


1 QO+0.04 


147.7/148 


RCS2318+0034 


mtl 


0.861H^ 


29li 


9 rf^+0.02 


304.8/300 


RCS2319+0030 


1131? 


n 54+0-03 

'-'•^^-0.02 


19.71^:? 


1 c«+0.03 


209.2/155 


RCS2319+0038 


lOOle' 


n fi5+0 03 


46+3 


9 «1+0.02 
^•"-■--0.02 


153.6/155 


RCS2320+0033 


1171^ 


n 61+0 03 
'-'■"-'^-0.02 


18.911:° 


1 74+0.03 

-■-■'^-0.03 


162.8/147 



'^Surface brightness / in units of 10 ^ photons sec ^ cm ^ arcsec ^ 
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Table 4. Integrated Spectral Fits (A = 2500) 



Cluster 




kT 


Z 




Nh 


yVDOF 

/V, / 




\h.Z^ kucl 


[keVl 


fZnl 




[10^° cm-2] 




RCS0224-0002 


329li 




0.3 




2.91 


48.7/65 






5 1+1-2 
^•-■--o.g 


0.21°: 


.4 
,2 




48.7/64 


RCS0439-2904 


123+?^ 


1 c;+0.3 
-■-•^-0.2 


0.3 




2.63 


8.0/10 


RCSl 107-0523 


29611 


4 2+0.8 


0.3 




4.24 


52.2/61 






4 2+0-6 

^■^-0.5 


0.71°: 


,5 
,3 




50.1/60 


RCS1326+2903 


2021^^ 


^■"J-o.s 


0.3 




1.16 


9.3/11 


{z = 1.01) 


1281^? 


-'-•'-'-0.3 


0.3 






8.9/11 


RCS1419+5326 


356l}3^ 


^•'-'-0.3 


0.3 




1.18 


126.3/125 






^ _i_n /I 

4 fi+o-4 
^•"-0.3 


0.31°: 


, 1 

.1 




126.2/124 


RCS1620+2929 


2701^^ 


Q q+1-3 
"-••^-o.g 


0.3 




2.72 


14.1/20 






S q+i-2 

"J-^-O.S 


0.31°: 


.8 
.3 




14.1/19 


RCS2318+0034 


4101^? 


O--'--0.9 


0.3 




4.13 


48.9/68 






5 8+1-2 


0.61°: 


,3 
,3 




47.4/67 


RCS2319+0030 


31911 


6.5+i 


0.3 




4.13 


44.8/34 






6li 


0.61°: 


.4 
.4 




44.3/33 


RCS2319+0038 


35111 


6.2l°| 


0.3 




4.16 


72.4/81 






K q+0.8 
'-'•^-0.7 


0.51°: 


,2 
,2 




70.9/80 


RCS2320+0033 


3231^^ 


5.9+2 


0.3 




4.14 


31.2/32 






6.01^ 


0.31°: 


,4 
.3 




31.1/31 



Note. — Single temperature fits within R25oo- When possible, a second 
fit was performed allowing both the temperature and abundance to vary. 
These fits are reported in the second fine (where there is one) for each 
cluster. In the case of RCS1326-I-2903, the second line indicates the result 
of fitting the integrated spectrum with a fixed abundance and a redshift 
of z = 1.01. 
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Table 5. Cluster Richness and Luminosity 



dluster 


-'-'gc,red 










[10^^ erg s-i] 


[10^^ erg s-i] 


RCS0224-0002 


945 ± 210 


9 1 +0.3 




RCS0439-2904 


1590 ± 460 


1 p:+0.5 




RCS1107-0523 


899 ± 280 


9 q+0.3 

Z.Og 2 


q c+0.3 


RCS1326+2903 


381 ± 275 


n 4+0-3 


-1 -1+0.5 
-'-•-'-0.5 


{z = 1.01) 


2670 ± 671 




9 7+1.0 
^- ' 1.0 


RCS1417+5305 


1879 ± 464 


1 q+0.9a 
-'-•'-'2.0 




RCS1419+5326 


1173 ± 224 






RCS1620+2929 


906 ± 236 


9 q+0.7 
^•'-'0.3 


q q+0.5 
O.Oq 7 


RCS2112-6326 


1011 ±400 


1 9+0.8a 
-■-■^2.5 




RCS2156-0448 


481 ± 166 


n 1 +0.1a b 
'-'■-'-0.3 




RCS2318+0034 


996 ±217 


U.Uq 4 


q+0.9 


RCS2319+0030 


1150 ±281 


O.Uq 4 


7 q+0.7 


RCS2319+0038 


1515 ±323 


' -"0.4 




RCS2320+0033 


578 ± 202 


4 9+0.5 


K q+0.5 
'-'■^0.6 



'^Bolometric X-ray Luminosity within 500 hjQ kpc, assuming 
a temperature of 4 keV 



^ID uncertain (IBarrientos et al. 1120071 ) 
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Table 6. Dynamical Comparisons 



Cluster 


cr 


]^g2.49 rjO.eSa 




[km s~"^] 


[km s""*"] 


RCS1107-0523 


700 ± 300 


785^?^ 



RCS1620+2929 
RCS2319+0038 



1050 ± 340 
860 ± 190 



7481}?^ 
10121^? 



Xue fc Will torn 
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Table 7. Mass Estimates (A = 2500) 



Cluster no Mgas M250O fgas 

[10-2 cm-3] [10^^ M0] [10^^ Mo] 



RCS0224-0002 
RCS0439-2904 
RCS1107-0523 
RCS1326+2903 

(z = 1.01) 
RCS1419+5326 
RCS1620+2929 
RCS2318+0034 
RCS2319+0030 
RCS2319+0038 
RCS2320+0033 



0.329 
0.545 



hO.009 
-0.010 
+0.016 
-0.016 
1 72+0.058 

n QOQ+O.OIO 
U.OZO_o 010 

0.4321H11 

A 07+0.047 
Z.^Z/_0 048 

U.O^O_oo35 

hO.Oll 
-0.012 
hO.019 
-0.019 
-1-0.033 
-0.033 
h0.019 
-0.019 



0.42- 



0.713: 
0.698: 
1.205 
0.699 



-0.05 
-0.05 

n7+o-oi 

U-U(_o.oi 



0.33 



-0.04 
-0.04 

o.iolHl 

0.061H1 
0.75lH^ 
0.291H? 

^■^ ' -0.06 

0.87t^:}^ 

1-0.07 
-0.07 



0.59 



12.151}:^^ 
0.86;°:}° 

7 SS+0.88 
' -^^-0.88 

07+0.54 

-0.56 

1 90-1-0.24 

11.391°:^^ 

-1-1.51 
1.45 
-2.82 
-2.92 
-1.34 
-1.35 
-1.59 
-1.60 
-2.06 
-2.10 



7.60 

24.2 
12.55 
16.54 
13.57 



0.0351°:°°^ 
0.0781°:°}J 

U.U^Z_o 005 

o.o34i°:r 

0.048inj^ 

U.UDO_o.oo5 

U.UO»_o.o08 

n 040+°-°°^ 

U.U^U_o.005 
U.UOZ_o Q06 

U.UOZ_o 006 

U-U++_o.o07 
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Table 8. Mass Estimates (A = 500) 



Cluster 


R-500 


Mgas 


M500 


fgas 




[kpc] 


[10^^ Mq] 


[10^^ Mq] 





RCS0224-0002 8191^°^ OAbSt^w 4.9751°;^^° 0.0311°;™^ 

RCS0439-2904 350^^^ 0.0471°;°°^ 0.3641°;°^^ 0.1291^?° 

RCS1107-0523 6651^^ 0.1131°;°}^ 1-7821°^^^ 0.0631°;°°^ 

RCS1326+2903 544+^}^ 0.0321°°°^ l-083t°-^-^^ 0.0301°°°^ 

(. = 1.01) UOtll' 0.0421°;°°^ 0.9121°;}^° 0.046^°;°^^ 

RCS1419+5326 802^^7 o.218in}9 2.600l°;}^^ 0.0841°;°°^ 

RCS1620+2929 627^^^ 0.09ll°:°l^ 2.3421°;^°^ 0.0391°;°°^ 

RCS2318+0034 979^^?=^ 0.199l°;°28 12.936;^:|^^^ 0.0151°;°°^ 

RCS2319+0030 749^^^^ 0.264l°:°3i s.igi+ojsg 0.083l°:°}2 

RCS2319+0038 8091^ 0.264l°:°| 4.007l°:|^ 0.0661^;^°^ 

RCS2320+0033 760^^^^ 0.184l°;°^5 2.8231^;^^^ 0.065l^:[^i^ 

RCS0439-2904 569t^^ O.lOOlg:^^^ 0.6191°;°?? 0.1611^;^^^ 



-52- 



Table 9. Fitting Parameters 









A 


= 2500 






A = 500 




Fit 




Sample 


Ci 




C'logY 


Ci 




ClogY 


E~'^Lx ■ 


-Tx 


RCS 


0.36 ±0.06 


2.05 ±0.34 


0.17 


0.59 ±0.05 


1.79 ±0.42 


0.15 








0.45 ±0.03 


2.0 (fixed) 


0.20 


0.65 ±0.03 


2.0 (fixed) 


0.19 






CNOC 


0.74 ±0.08 


2.31 ±0.31 


0.18 














0.85 ±0.01 


2.0 (fixed) 


0.19 












TOTAL 


0.56 ±0.07 


2.90 ±0.35 


0.28 














0.81 ± 0.01 


2.0 (fixed) 


0.33 








E-^Lx - 


EMtot 


RCS 


-0.03 ± 0.04 


1.38 ±0.12 


0.16 


-0.20 ±0.16 


1.03 ±0.28 


0.24 








0.06 ±0.02 


1.33 (fixed) 


0.19 


-0.28 ±0.02 


1.33 (fixed) 


0.33 






CNOC 


0.44 ±0.12 


1.26 ±0.21 


0.20 














0.48 ±0.01 


1.33 (fixed) 


0.23 












TOTAL 


0.07 ±0.10 


1.77±0.15 


0.29 














0.40 ±0.01 


1.33 (fixed) 


0.33 








EMtot - 


-Tx 


RCS 


0.29 ± 0.03 


1.48 ±0.27 


0.09 


0.76 ±0.08 


1.72 ±0.65 


0.22 








0.28 ±0.03 


1.5 (fixed) 


0.10 


0.69 ±0.03 


1.5 (fixed) 


0.23 






CNOC 


0.24 ±0.02 


1.83 ±0.13 


0.07 














0.27 ±0.01 


1.5 (fixed) 


0.08 












TOTAL 


0.28 ±0.02 


1.63 ± 0.18 


0.09 














0.27 ±0.01 


1.5 (fixed) 


0.08 








j-i— y/0 j- 


-Yx 


RCS 


0.32 ±0.05 


0.73 ±0.05 


0.11 


0.22 ±0.04 


0.65 ±0.10 


0.11 








0.40 ±0.03 


1.1 (fixed) 


0.22 


0.08 ±0.03 


1.1 (fixed) 


0.27 






CNOC 


0.50 ±0.08 


0.80 ±0.09 


0.12 














0.29 ±0.01 


1.1 (fixed) 


0.20 












TOTAL 


0.41 ±0.03 


0.88 ±0.04 


0.14 














0.30 ± 0.01 


1.1 (fixed) 


0.22 








eI"'m,,, 


-Yx 


RCS 


0.25 ± 0.03 


0.52 ± 0.05 


0.06 


0.43 ±0.08 


0.64 ±0.22 


0.23 








0.23 ±0.03 


0.581 (fixed) 


0.09 


0.37 ±0.03 


0.581 (fixed) 


0.25 






CNOC 


0.05 ±0.04 


0.63 ±0.05 


0.08 














0.04 ±0.02 


0.581 (fixed) 


0.09 












TOTAL 


0.18 ±0.02 


0.49 ± 0.03 


0.10 
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Table 9 — Continued 



A = 2500 A = 500 

Fit Sample Ci C2 ciogY Ci Ci criogY 

0.04 ±0.02 0.581 (fixed) 0.14 



Note. — Best fits to scaling relations cosmologically corrected by the 
factor E^. Temperature is in units of 5 keV; luminosity in 10^ erg s~^; 
mass in units of 10^^ Mq; in 4 x 10^^ M0 keV. Scatter along the 

1 /2 

Y-axis is calculated as [Ei=i,Ar (log - Ci - C2 log XJ^ /AT] ' . 
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Table 10. Fitting Comparisons 



Sample^ Ci C2 Redshift 







X — Tx 








A = 2500 












RCS 


0.36 ± 0.06 


2.05 ± 0.34 


0.6 < 


z 


< 1.0 


CNOC 


0.74 ± 0.08 


2.31 ±0.31 


0.1 < 


z 


< 0.6 


TOTAL 


0.56 ±0.07 


2.90 ±0.35 


0.1 < 


z 


< 1.0 


ASFOl 




2.08 ±0.06 


0.1 < , 


z 


< 0.45 


A = 500 












RCS 


0.59 ±0.05 


1.79 ±0.42 


0.6 < 


z 


< 1.0 


ETB04 


0.50 ±0.11 


3.72 ± 0.47 


0.4 < 


z 


< 1.3 




Ez-^Lx 


- EMtot 








A = 500 












RCS 


-0.20 ±0.16 


1.03 ±0.28 


0.6 < 


z 


< 1.0 


ETB04 


-0.63 ±0.32 


1.88 ±0.42 


0.4 < 


z 


< 1.3 



EMtot -Tx 



A = 2500 








RCS 


0.29 ±0.03 


1.48 ±0.27 


0.6 < ^ < 1.0 


CNOC 


0.24 ±0.02 


1.83 ±0.13 


0.1 < z < 0.6 


TOTAL 


0.28 ±0.02 


1.63 ±0.18 


0.1< z< 1.0 


APP05 


0.23 ± 0.05 


1.70 ±0.07 


z < 0.15 


ASFOl 


0.27 ±0.34 


1.51 ±0.27 


0.1<z< 0.45 


A = 500 








RCS 


0.76 ± 0.08 


1.72 ±0.65 


0.6< z< 1.0 


FRBOl 


0.52 ±0.45 


1.78 ±0.10 


z < 0.09 


APP05 


0.58 ±0.14 


1.71 ±0.09 


z < 0.15 


KV05 


0.51 ±0.31 


1.79 ±0.19 


0A< z < 0.7 


ETB04 


0.59 ±0.05 


1.98 ±0.3 


0.4 < ^ < 1.3 




Ez ^^^Lx-Yx 




A = 500 








RCS 


0.22 ±0.04 


0.65 ±0.10 


0.6< z < 1.0 


MOT 


-0.10 ±0.01 


1.1 ± 0.01 


0.1 < - < 1.3 


E'J'M,,, - Yx 



Table 10 — Continued 



Sample** 



Co 



Redshift 



A = 500 

RCS 
APP07 
KVN06 



0.43 ±0.08 
0.17 ±0.2 
0.27 ±0.006 



0.64 ±0.22 
0.55 ±0.03 
0.581 ±0.009 



0.6 < ^ < 1.0 
z < 0.15 
theory 



Re ferenced samples: (ASFOl nAUen. Schmidt. &: Fabian 



(2001) 
( 2005 ) 

(20071) 

mm 

(FRB Ol) 
EqOiI); 
(KVN06J 
(M07) 
fl2003l ^ 



(APP05) 
(APP07) 



Arnaud. Pointecouteau. fc Pratt 



Arnaud. Pointecouteau. fc Pratt 



(BMS04) 



Ettori et al. 



Borgani et al. . 
J2004ah: 



(ETB04) 

Finoguenov. Reiprich fc Bohringer 

( KV05) 



Kotov &: Vikhlinin 



Kravtsov. Vikhlinin fc ^ 



MaughanI fl2007f ): (SPF03) 



agai 




Sanderson et al. 
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Table 11. B geared Fitting Parameters 



Fit 


Sample 


Ci 


C2 


ClogY 


7" 




— iU.4U it z.oi 


o.Do it u.sy 


n Qo 






« 1 n -1- n n/i 

— D.iU It U.U4 


9 99 (f\vr^A\ 

z.zz (^nxcQj 


n 9/1 

U.Z4 






— 0.00 HI i.uu 


z.oi HI u.oo 


n 10 

u.oz 






K _i_ n no 
— O.DO It u.uz 


00 (f\vaA\ 

z.zz i^nxeQj 


n QO 

U.oZ 




TOTA T 


in 1 n -I- OQ 
— iu.iu It z.zy 


Q KS -I- n 7R 
O.Do It U. / D 


n /IQ 
u.4y 






^ -i- n no 

— O.DO It U.UZ 


Z.ZZ ^^IlxeQJ 


n /in 

U.4U 




1 CjUO 


— 4.4o It U. / 


1 S4 -1- n 04 

l.o4 It U.Z4 




rn 




— 0.04 it i.OD 


i.OD it U.04 


n 01 
U.zi 






_Q Qc; -u n DzL 
— 0.00 m u.u^ 


1.11 I iiXcU. I 


n 1 A 






— z.yu It u. ( 


1 nn -1- n 9/1 

i.uu It U.Z4 


n 1 /I 

U. 14 






ofi -1- n no 

— O.ZO It u.uz 


i . i i ^^nxeQ J 


n 1 

U. iO 




TOT A T 


/I ns -U n 0/1 
— 4.U0 it u.y4 


i.oo it u.oi 


n 1 fi 
u.io 






Q OS -1- n no 

— o.zo It u.uz 


111 (^v(^A\ 
1.1 ± ^^iixeuj 


n 1 fi 

U.ID 




1 rLiUo 


OQ -I- n /I 
— z.zy It U.4 


n 7s 4- n 1 Q 
U. / It U.io 




■1^^12500 




s on -1- Q 1 1 
— o.zU it o.ii 


9 fifi -I- 1 nR 

z.oD it i.UD 


n Qfi 
U.oD 






A OQ -L n n/i 
— 4.yo it U.U4 


1 <^7 /'f^-^^^^^ 
i.D/ (^nxeaj 


n 00 
U.zz 




CNOC 


-4.81 ± 1.51 


1.75 ±0.49 


0.26 






-4.66 ±0.02 


1.67 (fixed) 


0.27 




TOTAL 


-6.31 ± 1.36 


2.24 ±0.45 


0.30 






-4.74 ±0.02 


1.67 (fixed) 


0.28 


M200 


RCS 


-9.53 ±6.20 


3.54 ±2.11 


0.51 






-4.27 ±0.04 


1.67 (fixed) 


0.30 




CNOC 


-4.61 ± 1.53 


1.88 ±0.50 


0.26 






-4.06 ±0.02 


1.67 (fixed) 


0.26 




TOTAL 


-5.86 ± 1.43 


2.30 ±0.48 


0.32 






-4.12 ±0.02 


1.67 (fixed) 


0.28 




YE03 


-4.55 ±0.89 


1.64 ±0.28 






B07 


-5.70 ±3.4 


2.1 ± 1.2 
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Note. — Fits to richness scaling relationships. Lu- 
minosity is given in units of 10^^ erg s~^, tempera- 
ture in units of 5 keV, and mass in 10^^ Mq. Pa- 
rameters for the present work are measured within 
A = 2500 unless otherwise noted. Referenced 



sampl es are (ASFOl) 
)7) 

[E,=i,^ (log Y, - - C2 logX, 



Yee fc EUingsonl (120081 ) . and 



(B07) iBlindert et al. I (120071 ). Scatter (o"iogY) is given 

^2 



as 



/N] 



